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ABSTRACT 

We derive stellar metallicities, light-weighted ages and stellar masses for a magnitude- 
limited sample of 175,128 galaxies drawn from the Sloan Digital Sky Survey Data 
Release Two (SDSS DR2). We compute median-likelihood estimates of these param- 
eters using a large library of model spectra at medium-high resolution, covering a 
comprehensive range of star formation histories. The constraints we derive are set by 
the simultaneous fit of five spectral absorption features, which are well reproduced by 
our population synthesis models. By design, these constraints depend only weakly on 
the a/Fe element abundance ratio. Our sample includes galaxies of all types spanning 
the full range in star formation activity, from dormant early-type to actively star- 
forming galaxies. By analysing a subsample of 44,254 high-quality spectra, we show 
that, in the mean, galaxies follow a sequence of increasing stellar metallicity, age and 
stellar mass at increasing 4000A-break strength. For galaxies of intermediate mass, 
stronger Balmer absorption at fixed 4000A-break strength is associated with higher 
metallicity and younger age. We investigate how stellar metallicity and age depend on 
total galaxy stellar mass. Low-mass galaxies are typically young and metal-poor, mas- 
sive galaxies old and metal-rich, with a rapid transition between these regimes over 
the stellar mass range 3 x 10 9 < Af* <3x 1O 1O M . Both high- and low-concentration 
galaxies follow these relations, but there is a large dispersion in stellar metallicity 
at fixed stellar mass, especially for low-concentration galaxies of intermediate mass. 
Despite the large scatter, the relation between stellar metallicity and stellar mass 
is similar to the correlation between gas-phase oxygen abundance and stellar mass 
for star-forming galaxies. This is confirmed by the good correlation between stellar 
metallicity and gas-phase oxygen abundance for galaxies with both measures. The 
substantial range in stellar metallicity at fixed gas-phase oxygen abundance suggests 
that gas ejection and/or accretion are important factors in galactic chemical evolution. 

Key words: galaxies: formation, galaxies: evolution, galaxies: stellar content 



1 INTRODUCTION 

The ages and metallicities of stellar populations in nearby 
galaxies are direct tracers of the star formation and chemical 
enrichment histories of the Universe. Clues to the ages and 
metallicities of the stars may be inferred from the integrated 
spec tra of galaxies , using 'stellar population synthesis mod- 
els' fTinslevlll978l: IBruzual A.II1983: IBruzual A. fc Chariot) 
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iMarastorJ Il99l: IVazdekislll999l) . Analyses based on these 
models have been traditionally limited by the difficulty of 
deriving independent constraints on the age, star formation 
history, metallicity and dust content of a galaxy. For exam- 
ple, changes in age, metallicity and attenuation by dust all 
have similar effects on the colours and low-resolution spectra 
of galaxies. This gives rise to well-known 'near-degeneracies' 
in the constraints derived on these different parameters. 

The expectation is that these degeneracies may be 
broken, at least in part, by appealing to refined spec- 
tral diagnostics which are not sensitive to attenuation by 
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dust and have different sensitivities to age and metallic- 
ity. Studies in this area have focused on a set of 25 ab- 
sorption features defined and calibrated in the spectra of 
460 nearby Galactic stars obtained at Lick Observatory 
(e.g. lFabeill973l : IWorthev et al .11994 IWorthev fc Ottavianil 
1997). These studies all pertain to early-type galaxies, be- 
cause the lack of hot stars in the Lick library does not al- 
low the spectral interpretation of star-forming galaxies. The 
conclusion from these studies is that, for early- type galax- 
ies, comparisons of the strengths of metallic lines and age- 
sensitive Balmer lines can break the age-metallicity degener- 
acy, but only in a relative way: the derived ages and metal- 
licities appear to depend strongly on the specific choice of 
metal indices (e .g. iKun tschner et alJl200lt lEisenstein et alJ 
l2003l:lThomas et alfeOoUT This is because the Galactic stars 
used to calibrate the Lick indices have approximately solar 
metal abundance ratios at any metallicity, whereas the ra- 
tio of a-elements to iron is see n to increase from dw arf to 
massive early-type galaxies (e.g. IWorthev et al.lll992l) . 

A main weakness of the original calibration of Lick in- 
dices is that it relies on spectra which were not calibrated 
in flux, and for which the resolution (~ 9 A FWHM) is 
three times lower than achieved by modern spectroscopic 
galaxy surveys, s uch as the Sloan Digital Sky Survey (SDSS; 
lYork et aljEiooh . Thus, the high-quality spectra from these 
surveys must be degraded to the calibration and resolu- 
tion of the original Lick spectra for index-strength anal- 
yses to be performed. The situation has changed recently 
with the development of medium-high resolution (< 3 A 
FWHM), flux-calibrated population synth esis models in- 
cluding stars in the full t emperature range JVazdekisll200ll 
iBruzual fc Charlodl2003h . These models can be compared 
directly to high-quality observed spectra of both early-type 
and late-type galaxies. Several studies have also quantified 
the sensitivity of stellar absorption features to changes in 
element abundance ratios in galaxy spectra |rkntalo et alJ 
Il998l: iTrager et aO 1200(1 iThomas et al.ll2003bT) . These de- 
velopments coincide with the advent of large homogeneous 
samples of galaxy spectra gathered by modern surveys such 
as the SDSS. 

This is the first paper of a series in which we combine 
modern population synthesis techniques with the statistical 
power of the SDSS to investigate the connection between 
metallicity, age and stellar mass in nearby galaxies. Here, we 
u se the medium-high resolu tion population synthesis code 
of lBruzual fc Charlod feOOl hereafter BC03) to derive esti- 
mates of the metallicities, ages and stellar masses of a sample 
of ~ 2 x 10 5 nearby galaxies from the SDSS Data Release 
Two (DR2). We adopt a Bayesian statistical approach and 
derive full likelihood distributions for these physical param- 
eters by comparing the observed spectrum of each galaxy 
with a comprehensive library of model spectra correspond- 
ing to different star formation histories. The comparison is 
driven by the strengths of 5 spectral features selected to de- 
pend only weakly on the a/Fe ratio, which we measure in 
the same way in model and observed spectra. 

An important specificity of our work is that we de- 
rive the above constraints not only for quiescent, early- type 
galaxies, but also for late-type, star-forming galaxies, for 
which the contamination of stellar absorption features by 
nebular emission must be removed. We explore the relation- 
ships between metallicity, age and stellar mass and the de- 



pendence of these relationships on galaxy structure. In a 
companion paper, we exploit these results to investigate the 
physical origin of the colour-magnitude relation and of the 
relation between Mg2-index strength and velocity disper- 
sion for early-type galaxies. The total metal content of the 
local Universe and the distribution of metals as a function of 
galaxy properties will be the subject of a subsequent paper. 

The paper is organized as follows. In Section [5] below, 
we present our sample of SDSS spectra and the models used 
to interpret them. We also describe our method for deriving 
metallicities, ages and stellar masses from observed galaxy 
spectra. The results are presented in Section |HJ where we 
study the dependence of metallicity and age on total stellar 
mass and the age-metallicity relation as a function of galaxy 
structure. We also address in that section the influence of 
aperture bias on our derived parameters. We summarize and 
briefly discuss our results in Section 0] 



2 THE APPROACH 

In this section, we describe our approach for deriving esti- 
mates of light-weighted ages and metallicities from the ob- 
served spectra of SDSS galaxies. These spectra are discussed 
in Section 12.11 along with the models we use to interpret 
them. In Section [2.21 we select a set of spectral absorption 
features which we argue should most robustly constrain age 
and metallicity. Then, in Section ^. 31 we outline the statisti- 
cal approach we adopt to derive ages and metallicities from 
galaxy spectra. We also highlight the sensitivity of our con- 
straints on age and metallicity to the observational signal- 
to-noise ratio (S/N) of the spectra. 

2.1 Observed and model spectra 

The o bserved spectra we con sider are drawn from the SDSS 
DR2 jAbazaiian et alJl2004l . The SDSS is an imaging and 
spectroscopic survey of the high Galactic latitude sky, which 
will obtain u, g, r, i and z photometry of almost a quarter of 
the s ky and spectra of at least 700,000 objects iYork et alJ 
2000). The spectra are taken using 3"-diameter fibres, posi- 
tioned as close as possible to the centres of the target galax- 
ies. The flux- and wavelength-calibrated spectra cover the 
range from 3800 to 9200 A with a resolution of ~ 1800. 
Our sample, drawn from the SDSS DR2, includes unique 
spectra of 196,673 galaxies with Petrosian r-band magni- 
tudes in the range 14.5 < r < 17.77 (after correction for 
fo reground Galactic e xtinction using the extinction maps 
of lSchleeel et al.lfl998r> . The galaxies span the full range of 
types, from actively star-forming, late-type galaxies to dor- 
mant, early-type galaxies. The median redshift is 0.13. A 
more detailed description of this sample is given in Section 
13 below. 

To interpret these observed spectra in terms of physical 
parameters such as age, metallicity and stellar mass, we use 
the recent population synthesis models of BC03. These mod- 
els are based on 'STELIB', a newly avail able library of ob - 
served stellar spectra assembled bv lLe Borgne et alJ {2003). 
The models have a spectral resolution of 3 A FWHM across 
the whole wavelength range from 3200 to 9500 A. They are 
thus ideally matched to the SDSS spectra. 

We note that, to interpret the spectra of star-forming 



The Ages and Metallicities of Galaxies in the Local Universe 3 



and active galaxies with these models, we must first remove 
the contamination of the observed spectra by nebular emis- 
sio n lines. This is achiev ed using the procedure outlined 
bv lTremonti et all {2004), w hich is optimize d for use with 
SDSS galaxy spectra (see also lTremontl20o3) . This consists 
in performing first a non-negative least-squares fit of the 
emission-line-free regions of the observed spectrum, using a 
set of model template spectra broadened to the observed 
velocity dispersion (the template spectra correspond to 30 
instantaneous-burst models of different ages and metallici- 
ties computed using the BC03 code). Once the fitted spec- 
trum is subtracted from the observed spectrum, the residu- 
als can be fitted to Gaussian-broadened emission-line tem- 
plates. The method assumes a single broadening width for 
all the Balmer lines, and another (independent) width for 
all the forbidden lines. The strength of each line is fitted in- 
dependently. Then, the fitted emission lines are subtracted 
from the original observed spectrum to produce a 'pure' 
absorption- line spectrum suited to our analysis. 

2.2 Stellar absorption diagnostics of age and 
metallicity 

The strongest stellar absorption features in the opti- 
cal spectra of galaxies form the basis of the Lick sys- 
tem of spectral indice s <|Burstein et al.||l984|: |Gorgas et all 
[ 19931: IWorthev et all Il994l : IWorthev fe Ottavianil Il997t 
lTrageret^dT^)98t see also Section 0. Each index in this 
system is defined by a central 'feature bandpass' and two 
adjacent 'pseudo-continuum bandpasses'. The above studies 
have shown that some Lick indices are primarily sensitive to 
age, such as those based on H-Balmer lines, while others are 
primarily sensitive to metallicity, such as several Fe- and Mg- 
based indices at wavelengths between 4500 an d 5700 A. We 
also c onsider here the 4000- A break index of iBalogh et alJ 
( 1999), which we denote D4000, which is defined as the ra- 
tio of the average flux densities in the narrow bands 4000- 
4100 A and 3850-3950 A. 1 This index depends somewhat 
on metallicity but correlates more with the ratio of present 
to past-averaged star form ation rates in galaxies (see fig. 2 
of iKauffmann et al.ll2003al and fig. 27 of iBrinchmann et alJ 
120041) . 

In previous studies based on models with low spectral 
resolution, the strengths of Lick indices had to be mod- 
elled analytically as functions of the ef fective temperatures, 
gravities and metallicities of the stars I^ Vortjiev et al.lfl994t 
IWorthev fe Ottaviani|[l997t iGorgas et alJll999T) . These 'fit- 
ting functions' were not appropriate for hot stars, and hence, 
index-strength analyses had to be restricted to old stellar 
populations. Moreover, as mentioned in Section the stel- 
lar spectra on which the Lick indices were originally cali- 
brated had lower resolution than typical galaxy spectra to- 
day and were not flux-calibrated. These various weaknesses 
are resolved here by our adoption of the medium-high res- 
olution, flux-calibrated BC03 models. These models can be 
compared directly to the SDSS spectra of galaxies with any 
star formation history, and the Lick indices can be measured 
in the same way in model and observed spectra. 

Like previous models, however, the BC03 models rely on 

1 This index is sometimes denoted by D n (4000). 



a spectral library of nearby stars with near solar metal abun- 
dance ratios at any metallicity. The models are therefore 
expected to show discrepancies when compared to galax- 
ies where the abundance ratios differ from those of nearby 
stars. 2 In fact, such discrepancies appear to be responsible 
for the fact that some spectral features, such as CNi, CN2, 
TiOi, Ti02, Ca4227 and several Mg and Fe lines, were not 
well reproduced when BC03 compared their models with 
a sample of high-qu ality galaxy spectra dr awn from SDSS 
Early Data Release dStoughton et al.ll2002l) . 



Several studies have addressed the dependence of Lick 
index st rengths on chan ges in the relativ e rati os of heavy el- 
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Trager et all (2000: IVazdekisI l200ll: 


Thomas et al. 


2003b; 


Tantalo & Chiosill2004l: iThomas et all 


2004, see also 



BC03). These studies have led to the identification of com- 
posite Mg+Fe indices, which are sensitive to metallicity (i.e. 
the fraction by mass of all elements heavier than helium over 
the total gas mass) but show little sensitivity to a/Fe (i.e. 
the ratio of the total mass of a elements to the mass of iron; 
see footnote Among these, we use here 

[MgFe]' = V M S b (0-72 Fe5270 + 0.28 Fe5335), (1) 

as proposed bv lThomas et al.ll2003bl and 

[Mg 2 Fe] = 0.6 Mg 2 + 0.4 log(Fe4531 + Fe5015), (2) 

as defined in BC03. 

We wish to extract ages and metallicities from the SDSS 
galaxy spectra. We therefore fit simultaneously both metal- 
sensitive and age-sensitive indices. Among the 28 spectral 
indices studied by BC03, we concentrate on those that are 
best reproduced by the models. This requirement already 
excludes several Fe-based indices, the three Mg-based in- 
dices, Ca4227, Ca4455, the CN, TiO and NaD features (see 
fig. 18 of BC03). Among the remaining indices we identify 
those that are known to have at most a weak dependence 
on a/Fe. [Mg2Fe] and [MgFe]' are suitable metal-sensitive 
indices and H/3 is a suitable age-sensitive index. We also in- 
clude D4000, which is sensitive to the ratio of present to 
past-averaged star formation rate (see above). We are not 
aware of any study indicating the dependence of this index 
on a/Fe. Finally, to better constrain age we decided to in- 
clude also the two higher-order Balmer lines, H<5a and H7A. 
These indices have been recently shown to depend on varia- 
tions in eleme nt abundance ratios at metallicities aro und so- 
lar and above iThomas et alt 2004; Ko rn et alJ l2005L How- 
ever, we find that including them in our procedure does not 
produce results systematically different from those obtained 
without them, while it provides smaller errors on both age 
and metallicity estimates. Therefore, we decide to use also 
H<5a+H7a- We choose the sum of H<5a and H7A because 
it is better reproduced by the models than the two indices 
separately (see fig. 18 of BC03). Thus, our final set of in- 
dices is composed of D4000, H/3, RSa+R^a, [Mg 2 Fe] and 

2 The abundance ratio of a elements (such as N, O, Mg, Ca, 
Na, Ne, S, Si, Ti), which are produced mainly by Type II super- 
novae, to Fe-peak elements (such as Cr, Mn, Fe, Co, Ni, Cu, Zn), 
which are produced mainly b y Type la supernovae , is observed to 
vary in external galaxies (e.g., Worthe v et alll9 92: Thoma s et alJ 
l2003al) . 



4 A. Gallazzi et al. 



[MgFe]'. This is the minimum set that allows us to derive 
good constraints on metallicity and age simultaneously and 
to recover well the parameters of simulated galaxies (see sec- 
tion^] below). 

The solid histograms in Fig. show the distributions 
in 'resolving power' of these five spectral features for the 
galaxies in our sample. The resolving power is defined as 
the ratio between the 5-95 percent percentile range of the 
distribution of index strengths for all galaxies in the sam- 
ple, A/, and the observational error for each galaxy, crj. 
The resolving power is largest (~ 30) for D4000 and slightly 
lower (~ 10) for the other indices. Also shown as dotted his- 
tograms in Fig.Qare the analogous distributions in resolving 
power for a subsample of galaxies with mean signal-to-noise 
per pixel greater than 20. As expected, the distributions for 
these galaxies are shifted to higher values, because of the 
smaller observational errors. 



We note that the strengths of some spectral absorp- 
tion indices are sensitive to the ste l lar velocity dispersion 



in a galaxy (e.g. 


Davies et al.lll993t lLonehetti et al.lll998l 


iTraeer et alJll99£ 


; |Kuntschner||2004f). The indices that are 



most affected are those measured with the narrowest pseudo- 
continuum bandpass definitions. In particular this effect is 
seen in Fe-based indices and therefore also in the compos- 
ite Mg+Fe indices. The BC03 population synthesis code 
provides SSP spectra broadened to different velocity dis- 
persions. This allows us to compare each galaxy spectra 
with models that have a similar velocity dispersion (see Sec- 
tion 

2.3 Statistical estimates of age and metallicity 

We wish to estimate not only the most likely values of 
the ages and metallicities of SDSS galaxies, but also the 
accuracy of these values. To this purpose, we adopt a 
Ba yesian statistical approa ch, similar to the one outlined 
bv lKauffmann et al ] j2003aTl . The goal is to obtain the like- 
lihood distribution of a given parameter X in the space of all 
possible values of X. This is obtained by comparing the ob- 
servational data with a set of models that populate the space 
of all possible X according to a prior distribution which rep- 
resents our prejudice about the relative likelihood of differ- 
ent X values in absence of any data. 

We generate a library of 150,000 Monte Carlo realiza- 
tions of a full ran ge of physically plausible star formation his- 
tories. Following Kauffm ann et al] (1200381) . each star forma- 
tion history is parametrized in terms of a continuous model 
in which stars are formed from the time tf or m to the present 
according to the law oc exp(— ^t). Random bursts of 
star formation are superposed on this continuous model. We 
take the formation time tf or m to be uniformly distributed 
between 13.5 and 1.5 Gyr and the star formation time-scale 
parameter 7 to be uniform over the interval from to 1 
Gyr -1 . The bursts can occur at all times after tf OIUl with 
equal probability, set in such a way that 10 percent of the 
galaxies in the library experience a burst in the last 2 Gyr. 3 



They are parametrized in terms of the fraction A of stellar 
mass produced during the burst relative to the total mass 
formed by the continuous model. The ratio A is logarithmi- 
cally distributed between 0.03 and 4. During a burst, stars 
form at a constant rate for a time distributed uniformly 
in the range 3 x 10 7 -3 x 10 8 yr. The velocity dispersions 
of the models are distributed uniformly in the range 50- 
350 km s _1 . We further take the models to be distributed 
logarithmically in metallicity in the range 0.2-2.5 Zq and 
make the density of models drop smoothly as (log Z) 17/3 at 
metallicities from 0.2 down to 0.02 Zq, in order not to over- 
represent extremely metal-poor models. All stars in a given 
model have the same fixed metallicity, which we interpret as 
the '(optical) light- weighted' metallicity. 

For each model in the library we compute the following 
properties: 

(i) the strengths of the D4000, H/3, EW4+H7.4, [Mg 2 Fe] 
and [MgFe]' spectral indices, measured in the same way as 
in the SDSS spectra; 

(ii) the r-band light-weighted age, evaluated by the in- 
tegral t r = f* [dri/>(t - r) / r (r)r] / f* [dr^(t - r) f T (r)], 
where /, (r) is the total r-band flux produced by stars of 
age r. We refer below to the r-band light-weighted age sim- 
ply as the 'age' of a galaxy; 

(iii) the 2-band stellar mass-to- light ratio M*/L z , which 
accounts for the gas mass returned to the ISM by evolved 
stars (see also section 3.1 of BC03); 

(iv) the apparent u, g, r, i and z magnitudes of the model 
at redshifts between and 0.3 in steps of 0.01. 

The models in the library provide accurate simultane- 
ous fits to the strengths of the 5 spectral indices D4000, 
H/3, USa+H-Ja, [Mg 2 Fe] and [MgFe]' that we have selected 
to derive age and metallicity estimates from SDSS galaxy 
spectra. To illustrate this, Fig. shows the distribution of 
the differences between the best-fit and observed strengths of 
each index, in units of the observational error, for the galax- 
ies in our sample. Because of the dependence of the com- 
posite Mg+Fe indices on stellar velocity dispersion, when 
fitting the index strengths of an observed galaxy, we only 
include those models for which the stellar velocity disper- 
sion is within ±15 kms -1 of the observed one. 4 The plain 
histograms in Fig. [5] show the distributions for the sample 
as a whole, while the hatched histograms show the distribu- 
tions for those galaxies with a median S/N per pixel greater 
than 20. For both samples, the deviations between observa- 
tions and best-fit models are within the observational errors 
for all indices, as indicated by the comparison with a Gaus- 
sian distribution of unit standard deviation (dotted line). 
There is no strong correlation between the residuals of the 
different indices. 

The comparison of the strengths of D4000, H/3, 
H<5a+H7a, [Mg2Fe] and [MgFe]' in the spectrum of an ob- 
served SDSS galaxy with the strengths of these indices in ev- 
ery model spectrum in the library allows us to construct the 
probability density functions (PDFs) of physical parameters 



3 In [Kauffman n et al this parameter was set such that 

50 percent of the models in the library experienced a burst in 
the last 2 Gyr. We reduce this fraction to 10 percent because 
it provides a more uniform distribution of the models in light- 



weighted age. The influence on our results from changing this 
fraction is discussed in Section l2.4.2l 

4 The median uncertainty on velocity dispersion is ~ 15 kms -1 
for the full sample and only ~ 7 kms -1 for galaxies with S/N>20. 
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Figure 1. Distributions in 'resolving power' of the 5 spectral features selected to constrain the ages and metallicities of SDSS galaxies, 
as indicated. The resolving power is defined as the ratio between the 5%— 95% percentile range A/ of index strengths in the sample and 
the observational error aj of each galaxy. In each panel, the solid histogram shows the distribution for our sample of 196,673 SDSS-DR2 
galaxies, and the dotted histogram the distribution for the subsample of 44,347 galaxies with mean signal-to-noise per pixel greater than 
20. 
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Figure 2. Simultaneous fit of the 5 spectral features chosen to derive age and metallicity estimates from SDSS galaxy spectra using the 
model library described in Section 12.31 The histograms show the distributions of the differences between the best-fit (/ fit ) and observed 
(7 obs ) strengths of each index (as indicated), in units of the observational error (uj). The plain histograms are for the full sample of 
196,673 galaxies, while the hatched histograms are for the subsample of 44,347 galaxies with mean S/N per pixel greater than 20. Only 
those models for which the stellar velocity dispersion is within ±15 km s _1 of the observed one are included in the fit for each galaxy. 
For reference, the dotted line in each panel shows a Gaussian distribution with unit standard deviation. 



(such as age, metallicity and mass-to-light ratio) for that 
galaxy. This is achieved by assigning to each model a weight 
w = exp(— x 2 /2), where \ 2 ls calculated by comparing the 
strengths of the 5 indices measured in the observed spec- 
trum with those measured in the model spectrum, given the 
observational measurement errors. The PDF of a selected 
physical parameter is then simply given by the distribution 
in that parameter of the weights w of all the models in the 
library. 

The PDF of a given parameter can then be character- 
ized by its mode, which corresponds to the most likely value 
of the parameter, its median, which can differ from the mode 
for non-symmetric distributions, and by a confidence inter- 
val within which the parameter is constrained at a certain 



probability level. We often quote below the 68 percent con- 
fidence interval corresponding to the 16%-84% percentile 
range of the PDF, which would be equivalent to the ±lcr 
range for a Gaussian distribution. 

2.4 Accuracy of the estimates 

2.4-1 Dependence on observational properties 

We now want to illustrate the kind of constraints that can 
be obtained on the ages and metallicities of SDSS galaxies 
with d ifferent spectral propertie s using the method outlined 
above. iKauffmann et all <l2003al) have shown that the com- 
bination of H-Balmer lines and D4000 is a discriminating 
diagnostic of the recent star formation activity in galaxies. 
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We therefore draw 4 galaxies from our sample with high- 
quality spectra (median S/N per pixel larger than 30) lo- 
cated at different positions along the sequence occupied by 
SDSS galaxies in the HSa+Ht/a versus D4000 diagnostic di- 
agram. 

The solid distributions in Fig^show the constraints ob- 
tained on the metallicities (left-hand plot) and ages (right- 
hand plot) of these galaxies in two cases: when including 
only the age-sensitive indices D4000, H/3 and H<5a+H7a to 
constrain the fits (bottom panels), and when including also 
the metal-sensitive indices [Mg2Fe] and [MgFe]' (top pan- 
els). As expected, the ages are well constrained by D4000, 
H/3 and H5a+H7a alone, and the corresponding PDFs do 
not change appreciably when including also the constraints 
from metal-sensitive indices. In contrast, the metallicities 
are well constrained only when the metal-sensitive indices 
are included in the fit. In each panel in Fig. [3] the arrows 
indicate the median (longer one) and the 16 th and 84 th per- 
centiles (shorter ones) of the PDF. The 68 percent confi- 
dence interval becomes narrower when all the five indices 
are included, and the median of the distribution in this case 
is consistent with that obtained when including only the 
age-sensitive indices. We note that age and metallicity ap- 
pear to correlate with D4000 (indicated on the figure) for 
these 4 galaxies, and that the constraints on metallicity are 
weakest for the galaxy with the lowest D4000. This is not a 
coincidence, as we shall see in Section |21 below. 

It is important to note that the observed S/N directly 
influences the uncertainties in the age and metallicity esti- 
mates. To investigate this, we consider a set of galaxies with 
index strengths similar to those of the galaxies in Fig. |^1 but 
with lower median S/N per pixel, e.g. 20. We take the aver- 
age errors in the index strengths of these low-S/N galaxies 
to be typical of errors that would be measured at the same 
S/N for the galaxies in Fig. [3] We then mimic 100 real- 
izations of 'degraded spectra' of the high-S/N galaxies by 
randomly drawing index strengths from Gaussian distribu- 
tions centered on their index values, of widths the typical 
errors obtained from the lower-S/N galaxies. The PDFs of 
age and metallicity may be computed for each realization, 
and the average PDF of the 100 realizations is a good esti- 
mator of the results that would be obtained at low S/N for 
these galaxies. 

In the upper panels of Fig. [3] we show the average PDFs 
obtained in this way when the median S /N per pixel of each 
galaxy is degraded to 30 (long-dashed line), 20 (short-dashed 
line) and 10 (dotted line). A median S/N per pixel of 10 is 
not sufficient to constrain metallicity well (except for the 
most metal-rich galaxy, which has the spectrum with the 
strongest absorption features). As we increase the S/N, the 
distributions narrow down and converge to the PDFs ob- 
tained from the original high-S /N spectra. We conclude that 
a median S/N per pixel of at least 20 is required to reliably 
constrain metallicity. This is higher than the median value 
of ~ 15 for the SDSS-DR2 sample. The age estimates do not 
appear to be significantly affected by low S/N (except for 
the youngest galaxy) , probably because they are constrained 
strongly by D4000, which has small errors. The results of this 
exercise are summarized in Tabled 

The results of Fig. |3] would not change dramatically if 
we restricted the fit to only a subset of the selected spectral 
indices. We have checked this by using only H/3 to derive 



Table 1. Median-likelihood estimates of metallicity and age for 
the four galaxies shown in Fig.[!J] obtained by degrading the actual 
signal-to-noise ratio of the spectrum (S/N b s ) to 30, 20 and 10. 
The quoted errors are one half the 68 percent confidence range of 
the PDFs in the upper panels of Fig. EI Each column corresponds 
to a galaxy, ordered from left to right as in Fig. [3] 



log(Z/Z ) 


S/N = 


10 


-0.43±0.57 


-0.42±0.43 


-0.07±0.27 


0.09±0.12 


S/N = 


20 


-0.36±0.44 


-0.39±0.27 


-0.02±0.21 


0.12±0.08 


S/N = 


30 


-0.34±0.35 


-0.37±0.17 


0.00±0.16 


0.13±0.06 


S/N obs 


-0.41±0.22 


-0.38±0.12 


0.00±0.12 


0.14±0.04 


l°g(W»/r) 


S/N = 


10 


9.19±0.21 


9.68±0.18 


9.84±0.15 


9.92±0.11 


S/N = 


20 


9.12±0.15 


9.68±0.15 


9.83±0.14 


9.92±0.10 


S/N = 


30 


9.12±0.12 


9.66±0.13 


9.79±0.14 


9.92±0.09 


S/N obs 


9.12±0.08 


9.66±0.13 


9.79±0.13 


9.92±0.07 



age estimates for the galaxies in Fig. |3 and only [MgFe]' 
to derive metallicity estimates. The median-likelihood esti- 
mates of age and metallicity obtained in this way fall within 
the la errors of those shown in Fig. El The main effect of 
reducing the number of indices included in the fitting is to 
broaden the derived likelihood distributions. For the galax- 
ies of Fig. the average error on age increases from 0.10 
dex to 0.14 dex when fitting H/3 alone 5 , while the average 
error on metallicity increases from 0.13 dex to 0.29 dex when 
fitting [MgFe]' alone. 

2.4-2 Possible systematic uncertainties 

To identify any potential bias in our method, we also tested 
how well it can recover ages and metallicities of model galax- 
ies for which these parameters are known. We selected 3000 
models at random from the library described in Section ^. 31 
We added Gaussian noise to the index strengths of these 
models to reflect the average observational errors of our 
SDSS sample. The ages and metallicities recovered by our 
method for these models showed no systematic deviation 
from the true values. For both age and metallicity, the de- 
viations A log Z and A log t r for the 3000 models followed 
Gaussian distributions centred on zero of width ~ 0.15. In- 
terestingly, the deviations in age and metallicity appear to 
correlate with each other with a slope A log Z/ A log t r ~ 
—0.74. This is consistent with the age-metallicity degener- 
acy identified in Section 13.31 below. We have checked that 
the way in which A log Z and A log t r respond to (noise- 
induced) changes in metal-sensitive and age-sensitive indices 
for galaxies with different star formation histories cannot 
lead to spurious correlations between index strengths and 
metallicity /age residuals. 

We mentioned above that higher-order Balmer lines 

5 It has been noticed that higher-order Balmer lines can yield 
younger ages than those predicted by H/3, presum ably because 
of the dependence of h igher-order lines on a/Fc (e.g. Kuntschncr 
l2000t lEisenstein et al .112003?) . However, we do not see any signif- 
icant difference in the results when using either H<5a+H7a only 
or H/3 only. 
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Figure 3. Probability density functions of stellar metallicity (left-hand plot) and r-band light-weighted age (right-hand plot) for 4 SDSS 
galaxies with high-quality spectra (median S/N per pixel larger than 30) and different D4000 and H<5a+H7a strengths (indicated in the 
bottom panels of the left-hand plot). The solid PDFs in the bottom panels were obtained when including only the age-sensitive indices 
D4000, H/3 and H(5a+H7a to constrain the fits. Those in the top panels were obtained after including also the metal-sensitive indices 
[Mg2Fe] and [MgFe]'. In each panel, the arrows indicate the median (longer one) and the 16 th and 84 th percentiles (shorter ones) of 
the PDF. The long-dashed, short-dashed and dotted PDFs in the top panels show the constraints obtained when degrading the original 
galaxy spectra to a median S/N per pixel of 30, 20 and 10, respectively (see text for detail). 



may be sensitive to the a/Fe ratio. To quantify the potential 
error on our metallicity and age estimates for galaxies with 
enhanced a/Fe relative to solar , we compared the predic- 
tions of the lThomas et all i2004) stellar populations models 
at solar element abundance ratio, a/Fe = 0, with those for 
a/Fe — 0.3 (a typical ratio for massive elliptical galaxies). 
We chose 3000 models at random from our library and per- 
turbed their index strengths with Gaussian noise, as above. 
In addition, we increased the values of FLSa and H7A to re- 
flect the difference between the a/Fe = and a/Fe = 0.3 
tracks, interpolating in metallicity and age to the values 
for each model considered. The distributions in A log Z and 
Alogt r are well represented by Gaussian centred on +0.05 
and —0.05, respectively, of width given by the average error 
on metallicity and age (0.2 and 0.13 respectively). This test 
indicates that we tend to overestimate the stellar metallic- 
ities and underestimate the light-weighted ages of galaxies 
with supersolar abundance ratios by ~ 0.05 dex. However, 
this seems in contradiction with the fact that the results 
on the galaxies in our sample do not vary systematically 
if we include or exclude H<5a+H7a in the fit. This is true 
also for massive early-type galaxies, w hich are likely to be 
a-enhanced (e.g. IWorthev et aljfl99a) . The offset of ~ 0.05 
dex may thus be regarded as an upper limit to the error in 
the ages and metallicities derived from our analysis in the 
case of non-solar abundance ratios. 

Another possible source of systematic error is the choice 
of prior according to which our model library populates the 
parameter space. In particular, the mix of continuous and 
bursty star formation histories may influence the physical 
parameters in which we are interested, mainly the light- 
weighted age. To test for this effect, we generated a Monte 



Carlo library with a modified prior, by increasing to 50 per- 
cent (instead of 10 percent) the fraction of models that can 
undergo a burst of star formation in the last 2 Gyrs. We 
then compared the (median-likelihood) estimates of light- 
weighted age and stellar metallicity derived with this mod- 
ified prior and with our standard prior. Increasing the frac- 
tion of bursts, we derive ages ~0.07 dex younger, on aver- 
age, than those derived with our standard prior. This bias 
mainly affects old, early-type galaxies. Similarly the metal- 
licities are on average 0.04 dex higher than those derived 
with the standard prior. The effect of these offsets on the 
relations between age, metallicity and stellar mass that we 
discuss in Section 13.21 is very small. The zeropoint of the 
relations changes according to the offsets in metallicity and 
age reported above, but the shape of the relations remains 
identical. 



2.4-3 Stellar mass estimates 

In the remainder of this paper, we will be interested in the 
dependence of age and metallicity not only on directly ob- 
served properties, such as spectral features and morphol- 
ogy, but also on ' derived ' quantities, such as stellar mass. 
Kauf fmann et all (l2003al) derived Bayesian likelihood esti- 
mates of the stellar masses of a sample of SDSS galaxies, 
based on fits of the H5a and D4000 absorption indices. We 
use a similar approach here and estimate stellar masses for 
all the galaxies in our sample, based on the fits of D4000, 
H/3, RSa+FL-ja, [Mg 2 Fe] and [MgFe]'. We compute the PDF 
of the stellar mass for each galaxy by scaling the 2-band 
mass-to- light ratio M*/L z of each model to the observed, 
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total z-band luminosity of the galaxy. 6 We compute M*/L z 
at the observed galaxy redshift and include the effects of 
attenuation by dust. The z-band attenuation A z is inferred 
from the difference between the emission-line corrected fi- 
bre r — i colour 7 of the galaxy and the r — i colour of the 
redshifted (dust-free) model, assuming a single p ower law 
(cx A -0,7 ) attenuation curve JCharlot fc Falll200ch . 

The true attenuation can of course not be negative, 
but imposing strictly A z > would not account for the er- 
rors affecting the measurements of indices and magnitudes. 
Moreover, at high metallicities, there is a potential prob- 
lem that the dust-free models providing the best fits to the 
observed absorption-line strengths can be redder than the 
observed galaxy (implying negative apparent attenuations). 
This problem arises because of a discrepancy between the 
spectral absorption features and the integrated colours of 
the BC03 models at the metallicity 2.5 Zq. The colours of 
these models correspond to stellar evolutionary tracks with 
(and colour-temperature calibrations for) Z = 2.5 Zq, but 
the absorption- line strengths had to be calibrated using stel- 
lar spectra for slightly lower metallicity, 1.6 < Z/Zq < 2 (see 
Appendix A of BC03). Hence, the metallicity scale, which is 
linked to the evolutionary tracks and colours, may be biased 
high at supra-solar metallicities. This problem affects mainly 
the most metal-rich, early-type galaxies. To account for this 
bias, we include models producing dust attenuations down 
to A z — —0.1 when computing the PDFs of stellar mass, 
stellar metallicity and age. 8 

A com parison of our ste l lar ma ss estimates with those 
derived bv lKauffmann et alJ i2003al which were based on 
a preliminary version of the BC03 models and an early cal- 
ibration of the SDSS spectra) for the galaxies in common 
between the two samples shows overall consistency, with a 
scatter of ~ 0.16 dex. Our stellar mass estimates tend to be 
systematically larger than theirs by ~ 0.1 dex. A difference 
of almost +0.04 dex can be attributed to the different prior 
used to generate the Monte Carlo library. The remaining 
offset is likely to originate from their exclusion of all models 
with r — i colour redder than that observed, while we include 
models which imply A z 'attenuations' down to —0.1. 



3 THE AGES AND METALLICITIES OF 
NEARBY GALAXIES 

We now use the models described in Section [5] to estimate 
ages and metallicities f or a sample of 196 ,673 galaxies drawn 
from the SDSS-DR2 dAbazaiian et alJ l20o3 . see also Sec- 
tion above). Fig. 2] shows the distributions in median 

6 We assume that the M*/L z ratio is the same for the whole 
galaxy as it is in the region sampled by the fibre. See Fig. H3l for 
the typical fraction of light that enters the fibre. 

7 Fibre magnitudes are obtained using the SDSS photometry di- 
rectly out to the radius covered by the fibre. They are thus directly 
comparable to quantites derived from spectroscopy, provided that 
the fibre is positioned at the centre of the galaxy. Corrections for 
emission lines are obtained by comparing the magnitudes mea- 
sured off the spectrum before and after removing emission lines. 

8 The difference in r—% colour between two old stellar populations 
of metallicities 2.5 and 2Zq, when interpreted as a colour 
excess, corresponds to a z-band attenuation A z ~ 0.1. 



S /N per pixel (upper panel) and in redshift (lower panel) of 
this sample. Although the distribution in redshift extends 
up to 0.3, in what follows we discuss results only for galax- 
ies in the redshift range 0.005 < z 0.22. The lower limit 
is intended to avoid redshifts for which deviations from the 
Hubble flow can be substantial, but still allows us to include 
galaxies at very low luminosity (corresponding roughly to a 
lower mass limit of ~ 1O 8 M0). The upper limit corresponds 
roughly to the redshift at which a typical 10 1 Mq galaxy is 
detected with median S/N per pixel greater than 20. These 
cuts leave us with 175,128 galaxies. 

In some applications below, it will be useful to distin- 
guish between different morphological types of galaxies. This 
can be achieved on the basis of the 'concentration param- 
eter' C = R90/R50, defined as the ratio of the radii en- 
closing 90 and 50 percent of the Petr osian r-band luminos- 
ity of a galaxy. IStrateva et alJ fcOOll) and IShimasaku et all 
(2001) have shown that the concentration parameter al- 
lows a rough c lassification into galaxy morphological types. 
IStrateva et alJ i200ll) propose a cut at C = 2.6 to sepa- 
rate early- from late-type galaxies. To limit the contamina- 
tion between the two types, here we define those galaxies 
with C ^ 2.8 to be 'early-type' and those galaxies with 
C ^ 2.4 to be 'late-type'. The dotted and dashed lines in 
Fig- El show separately the distributions in median S/N per 
pixel and in redshift of the resulting subsamples of early- 
and late-type galaxies, respectively. Early-type galaxies gen- 
erally have a higher median S /N per pixel and are detected 
out to higher redshifts than late-type galaxies, as expected 
from their higher surface brightnesses. 

We have shown in Section 12.4.11 and Fig. |H] the im- 
pact of the signal-to-noise ratio on our ability to derive 
reliable constraints on age and, especially, metallicity. In 
the remainder of this paper we focus on the subsample of 
44,254 galaxies with median S/N per pixel greater than 20, 
for which the constraints on metallicity are the most ac- 
curate. The properties of these high-S/N galaxies are sum- 
marized in Fig. |K| where we show their distribution (shaded 
histogram) in z-band surface brightness (averaged within 
the Petrosian R50 radius), Petrosian r-band absolute mag- 
nitude, redshift and D4000, compared to the distributions 
for the original sample (dot-dashed line). The red and blue 
lines distinguish early- type from late- type (high-S/N) galax- 
ies. The S/N requirement excludes about 75 percent of our 
sample. As can be seen from Fig. 03 this biases the sample 
toward higher-surface brightness, more concentrated, lower- 
redshift galaxies. There is not a bias instead in the lumi- 
nosity distribution. In the redshift range 0.005 < z ^ 0.22, 
high-S/N galaxies typically have z-band surface brightnesses 
brighter than 21 mag arcsec -2 . They lie primarily at red- 
shifts less than 0.15, where they account for about 40 per- 
cent of the galaxies. They also represent roughly 40 percent 
of the high-concentration sample but only 10 percent of the 
low-concentration sample. Similarly, they constitute about 
30 percent of the galaxies with D4000 > 1.6 and 15 percent 
of those with D4000 < 1.6. Our selection in S/N therefore 
excludes a substantial fraction of diffuse systems, with po- 
tentially subsolar metallicities. We will turn our attention to 
these systems in the third paper of this series. Their omis- 
sion here does not alter significantly the results presented 
in the next sections. In particular, we have checked that the 
median ages and metallicities found for galaxies with stellar 
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Figure 4. Distributions in median S/N per pixel (upper panel) 
and in redshift (bottom panel) of the sample of 196,673 SDSS- 
DR2 galaxies (solid histograms). In each panel, the dotted and 
dashed histograms show the corresponding distributions for the 
subsamples of high-concentration (C 2.8), early-type galaxies 
and low-concentration (C ^ 2.4), late-type galaxies, respectively. 



masses less than ~ 2 x 10 9 Mq in Sections 13. 2l and l3.3l below 
would change by at most 0.2 dex if we included galaxies with 
low S/N in our sample. 



3.1 Age, metallicity and mass distributions 

We present here the distributions in age, metallicity and 
stellar mass of the 44,254 galaxies in our final sample. We 
recall that we estimate r-band light-weighted ages but as- 
sume that all stars in a given galaxy have a single metal- 
licity (interpreted as the optical light-weighted metallicity). 
The left-hand panels of Fig. |S| show, from top to bottom, 
the distributions of the median-likelihood estimates of stel- 
lar metallicity, age and stellar mass for the galaxies in our 
sample (solid histograms). The distribution in metallicity 
peaks around \og(Z / Zq) ~ 0.1, the number of galaxies with 
sub-solar metallicities decreasing smoothly down to 10 per- 
cent of solar. The distribution in age indicates that most of 
the galaxies have fairly old mean stellar populations, with 
t r ~ 7 — 8 Gyr. The distribution in stellar mass reveals that 
our sample is dominated by relatively massive galaxies, with 
M* ~ 1O 11 M0, while only a small fraction of the sample has 
masses below 1O 1O M . 

The right-hand panels of Fig. |5] show the distributions 
of one half the 68 percent confidence range in the estimates 
of log(^/Z ), log(t r /yr) and log(M*/M ). Stellar mass is 
the best constrained parameter, with a typical random un- 
certainty of ~ 0.08 dex. Almost all the galaxies in the sample 
have a stellar mass estimate with an uncertainty less than 
0.1 dex. Also, the ages are constrained within 0.2 dex for 
almost all the galaxies, the typical uncertainty being ~ 0.12 
dex. Similarly the average error on stellar metallicity is 0.12 



dex, but the distribution extends to larger errors (~ 0.3) 
than for light-weighted age. 

As a comparison, the dotted histograms show the dis- 
tributions in the derived parameters and the associated un- 
certainties for the original sample of 175,128 galaxies. There 
is a significant tail of galaxies with very broad stellar metal- 
licity PDFs, resulting in uncertainties greater than 0.25 dex. 
This further demonstrates that the uncertainties in metal- 
licity estimates depend sensitively on the S/N in the ob- 
served spectra, as discussed in Section f2. 41 Also, the distri- 
bution in stellar metallicity shows that with our S/N cut we 
have excluded preferentially galaxies with subsolar metal- 
licity. High-metallicity galaxies tend to be associated with 
small errors and low-metallicity galaxies with large errors 
(presumably because of their weaker absorption lines). In- 
cluding low-S/N galaxies also increases the average error on 
light-weighted age to ~ 0.15 dex, but the trend of increasing 
error with decreasing age is less strong than the analogous 
trend for metallicity. 



3.2 Relations between age, stellar metallicity, 
stellar mass and gas-phase metallicity 

It is interesting to examine how age, metallicity and stel- 
lar mass (and the errors on these quantities) are distributed 
among galaxie s with different star formation histories. As 
shown bv lKauffmann et al ] J2003a|) . galaxies with different 
star formation histories populate different areas of a dia- 
gram defined by D4000 and the strength of a H-Balmer ab- 
sorption line (they illustrated this result for H<5a)- Galaxies 
with smooth star formation histories form a sequence ex- 
tending from actively star-forming galaxies (corresponding 
to small D4000 and strong H-Balmer absorption) to quies- 
cent early-type galaxies (corresponding to large D4000 and 
weak H-Balmer absorption) . Galaxies which experienced re- 
cent bursts of star formation exhibit the strongest H-Balmer 
absorption at fixed D4000. 

We now explore how the 44,254 high-S/N galaxies in 
our sample populate such a diagram as a function of their 
physical parameters derived in Section 13.11 above. We con- 
sider here for consistency the diagram defined by D4000 
and the HSa+'R^a absorption index. In the left-hand pan- 
els of Fig. |7] we have binned and colour-coded this diagram 
in order to reflect, from top to bottom, the average stellar 
metallicity, the average age and the average stellar mass of 
the galaxies falling into each bin. The widths of the bins 
in D4000 and USa+Ht/a correspond roughly to the mean 
observational errors in the two quantities for galaxies with 
median S/N per pixel greater than 20 (0.04 and 0.4, respec- 
tively) . 

The bottom-left panel of Fig.Qshows that stellar mass 
increases with D4000 along the sequence populated by SDSS 
galaxies in the H<5a+H7a versus D4000 diagram. The aver- 
age mass of galaxies with D4000 < 1.2 is less than 1O 1O M , 
while that of galaxies with D4000 > 1.9 reaches 10 n M Q . 
There is a smooth transition between these two regimes 
in the central region of the diagram. This confirms the 
trends in D4000 and H£a as funct ions of stellar mass ob- 
served bv Hvauffma nn et alj j2003bj) (see their fig. 1). Both 
stellar metallicity and age also appear to increase with in- 
creasing D4000 (upper- and middle- left panels of Fig. [7J. 
The variation in age is very smooth. In the intermediate 
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Figure 5. From top-left to bottom-right: distribution in z-band surface brightness, r-band absolute Petrosian magnitude (k-correctcd at 
z = 0), redshift and 4000A- break strength for the galaxies in our sample. The dot-dashed line represents the distribution for the original 
sample of 175,128 galaxies in the redshift range 0.005 < z sC 0.22. The shaded histogram describes the distribution for the final sample of 
44,254 high-S/N galaxies, while the red and blue histograms separate (high-S/N) early-type (C J? 2.8) and late-type (C $C 2.4) galaxies, 
respectively. (Note that the dot-dashed and the shaded histograms are both normalized to unit area.) 



regime, age appears to decrease with increasing TLSa+R^a 
at almost fixed D4000. The variation in metallicity is more 
noisy, partly because of the larger errors associated with 
metallicity estimates, especially for low D4000 (see Sec- 
tion 123- Near the middle of the galaxy sequence, i.e., for 
1.4 < D4000 < 1.8, stronger H-Balmer absorption appears 
to be associated to not only younger but also more metal- 
rich stellar populations. This is consistent with the idea 
that these galaxies could have experienced a burst of metal- 
enriched star formation about 1-2 Gyr ago. 9 

We have checked that the above trends in stellar metal- 
licity, age and stellar mass do not change significantly if we 
include galaxies with lower S/N values. The main effect is to 
increase the scatter in these physical parameters along the 
H<5a+H7a versus D4000 sequence, mainly because of the 
larger observational errors in these two indices. Moreover, 
the average S/N for the full sample does not vary signifi- 
cantly along the relation. For these reasons, we believe that 
our results are not strongly biased by our cut in S/N. 

In the right-hand panels of Fig.0 we have colour-coded 
the H5a+H7a versus D4000 diagram to reflect the average 
uncertainties in the determinations of stellar metallicity, age 
and stellar mass (from top to bottom) for the galaxies falling 
into each bin. These diagrams illustrate how the tightness of 



9 At fixed D4000 ~ 1.5, this trend in age and metallicity is not 
consistent with (and more pronounced than) the age-metallicity 
degeneracy described in Section l3.3l 



the constraints derived on the various physical parameters 
depends on the strengths of H<5a+H7a and D4000. The er- 
rors in all three parameters tend to be larger at lower D4000, 
the trend being especially strong for stellar metallicity. We 
emphasize that the right-hand panels of Fig. |7| do not re- 
flect the rms scatter in the various physical parameters, but 
rather the average uncertainties associated with determina- 
tions of these parameters. For reference, the rms scatter is 
typically comparable to (or smaller than) the average uncer- 
tainty for stellar metallicity and age, but always larger than 
the average uncertainty for stellar mass. 

The results of Fig. |7| suggest that both stellar metal- 
licity and age correlate with stellar mass. We investigate 
this further in Fig. |H] where we show the distributions of 
metallicity (panel a) and age (panel b) as function of stel- 
lar mass for the 44,254 high-S/N galaxies in our sample. 
Rather than assigning each galaxy its median-likelihood es- 
timate of each of the three parameters, we keep here the 
whole information contained in the PDFs. In Fig. |HJt,, the 
likelihood distribution of stellar metallicity as a function of 
stellar mass was obtained by coadding the normalized 2-D 
likelihood distributions of metallicity and stellar mass for 
all the galaxies and then re-normalizing along the metallic- 
ity axis in bins of stellar mass (we adopted a bin width of 
0.2 dex, roughly comparable to the 68 percent confidence 
range in stellar mass estimates) . The solid line indicates the 
median of the final conditional distribution, and the dashed 
lines the 16 th and 84 th percentiles. The percentiles of the 
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Figure 6. Distributions of the median-likelihood estimates of stellar metallicity (top-left panel), age (middle- left panel) and stellar mass 
(bottom-left panel) for 44,254 SDSS-DR2 galaxies with redshifts in the range 0.005 < z sC 0.22 and median S/N per pixel greater than 20 
(solid histograms). The right-hand panels show the corresponding distributions of errors, computed as one half the 68 percent confidence 
ranges in the estimates of log(Z/Z ), log(t r /yr) and log(M* /M ). The dotted histograms are for the original sample of 175,128 galaxies 
(all histograms are normalized to unit area). 



metallicity and age distributions at fixed stellar mass are 
provided in Table [5] 

The median shows that metallicity increases with stel- 
lar mass, from roughly 20 percent of solar for galaxies with 
masses below 1O 1O M to about 1.4 times solar for those 
with masses above 1O 11 M0. In between, stellar metallic- 
ity increases rapidly with stellar mass, the trend becom- 
ing shallower at M« > 1O 1O,5 M . Fig. [Ht> shows the anal- 
ogous conditional distribution of age as a function of stel- 
lar mass. The median (r-band light-weighted) age increases 
from roughly 1 Gyr for galaxies with masses around lO 9 A/ 
to about 6.5 Gyr for those with masses around 1O 1O,5 M0. 
For more massive galaxies, the median age continues to in- 
crease, but only by about 0.2 dex (~ 1.6 Gyr) over one order 
of magnitude in stellar mass. 

Although the median metallicity and age both increase 
with increasing stellar mass, the probability levels (indicated 
by the colour) and the 68 percent confidence ranges reveal 
broad distributions in the estimates of both parameters. We 
now explore whether the larger scatter in both distributions 
at smaller masses can be accounted for by the larger uncer- 
tainties in metallicity and age estimates for more metal-poor 
and younger galaxies, or whether it is indicative of an intrin- 
sically broad distribution in metallicity and age for low-mass 
galaxies. In the bottom panels of Fig. |H| we compare the 
scatter in these relations with the uncertainties in the es- 
timates of metallicity (panel c) and age (panel d) at fixed 
stellar mass. In each panel, the diamonds show the differ- 
ence between the 16 th and 84 th percentiles of the conditional 



Table 2. Median (P50) and percentiles (P16, P84) of the distri- 
butions in stellar metallicity and age as a function of stellar mass 
(Fig.&,b). 



log(Z/Z ) log(Wyr) 



log(Af./Af ) 


P50 


P16 


P84 


P50 


P16 


P84 


8.91 


-0.60 


-1.11 


-0.00 


9.06 


8.80 


9.46 


9.11 
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distribution shown above as a function of stellar mass. The 
stars show the mean 68 percent confidence range in the indi- 
vidual estimates of metallicity or age for the galaxies falling 
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Figure 7. Physical parameters of 44,254 SDSS-DR2 galaxies with median S/N per pixel greater than 20 as a function of position in the 
IWa+I^a versus D4000 diagram. In the left-hand panels, the diagram has been binned and colour-coded in order to reflect the average 
stellar metallicity (top), age (middle) and stellar mass (bottom) of the galaxies falling into each bin. In the right-hand panels, the diagram 
has been binned and colour-coded to reflect the average uncertainty (corresponding to one half the 68 percent confidence interval) in the 
derived stellar metallicities (top), ages (middle) and stellar masses (bottom) of the galaxies falling into each bin. The widths of the bins 
in D4000 and H^A+UfA correspond roughly to the mean observational errors in these quantities (0.04 and 0.4, respectively). 



into each stellar-mass bin. 10 Fig.|HJ; shows that the scatter in 
the mass-metallicity relation is always larger than the uncer- 
tainties in the metallicity estimates, by at least 0.1 dex, and 
is largest (compared to the error) at about 1O 1O M0. Simi- 
larly, the scatter in the mass-age relation (Fig. |HJi) is always 
larger than the uncertainties in the age estimates and is also 
largest at about 1O 1O M0. Fig. |H| clearly indicates that more 
massive galaxies are older and more metal-rich, while less 
massive galaxies are younger and more metal-poor. Metal- 
licity and age, however, are not uniquely determined by stel- 
lar mass. There is an intrinsic scatter in both parameters, 



10 The stars in Fig. |SJ;,d are not computed using the joint likeli- 
hood distribution, but rather by adopting the median of the M„ 
PDF as our estimate of stellar mass. The comparison of the stars 
with the diamonds is consistent, given that the width of the stellar 
mass bins is comparable with the error in stellar mass estimates. 



which is particularly evident at intermediate stellar masses. 
The results of Section r3.3l below indicate that part of this in- 
trinsic scatter can be accounted for by differences in galaxy 
morphology, the intermediate-mass regime corresponding to 
the transition between disc-dominated and bulge-dominated 
galaxies. However, a significant scatter persists even when 
considering the two classes of galaxies separately. 

iTremonti et al.l i2004) have shown that the emission- 
line galaxies in the DR2 sample exhibit a tight relation 
between stellar mass and gas-phase oxygen abundance, as 
determined from their nebular spectra. The gas-phase oxy- 
gen abundance increases steadily from low to high stel- 
lar masses, and then gradually flattens around masses of 
10 10 5 M W . Ther e is a s triking similarity between the relation 
ITremonti et al.l (120041) found for star forming galaxies and 
the relation found in Fig.[S]above between stellar metallicity 
and stellar mass, which includes both star-forming and qui- 
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Figure 8. Panels a and b: conditional distribution of stellar metallicity and age as a function of stellar mass for 44,254 SDSS-DR2 galaxies 
with median S/N per pixel greater than 20. Each distribution was obtained by coadding the normalized 2-D likelihood distributions of 
the desired parameter (stellar metallicity or age) and stellar mass for all the galaxies and then re-normalizing along the y-axis in bins of 
stellar mass (of width 0.2 dex). The solid line indicates the median of the final conditional distribution, and the dashed lines the 16 th 
and 84 th percentiles. Panels c and d: comparison of the scatter in the metallicity and age distributions with the mean uncertainty of the 
metallicity and age estimates as a function of stellar mass. The diamonds show the difference in \og(Z / Zq) or log(t r /yr) between the 
16 th and 84 th percentiles of the conditional distributions of panels a and b, while the stars represent the mean 68 percent confidence 
range in \og(Z / Zq) or log(t r /yr) for the galaxies in each stellar mass bin. 



escent galaxies. However, the dispersion in the stellar metal- 
licity versus stellar mass relation is much larger than that in 
the relation between gas-phase oxygen abundance and stel- 
lar mass (±0.1 dex). This may partly reflect the larger errors 
associated with stellar metallicity determinations. 

In Fig. |U] we plot stellar metallicity against gas-phase 
oxygen abundance, 12+log(0/H), for the 7462 high-S/N 
galaxies of our sample for which both measures are avail- 
able. We note that gas-phase metallicity is best determined 
for star-forming galaxies, whereas stellar metallicity is best 
determined for early- type galaxies. The small dots show 
the median-likelihood estimates of the stellar and gas-phase 
metallicities for each galaxy. The larger points indicate the 
median stellar metallicity in bins of 12+log(0/H), while the 
outer solid lines show the corresponding 16 th and 84 th per- 
centiles. There is a relation between the two parameters 
with approximately unit slope. However, the stellar metal- 



licity is generally lower than the gas-phase metallicity (by 
~ 0.5dex), as demonstrated by the one-to-one relation (dot- 
dashed line). This is expected since the nebular metallicity 
traces the metallicity of the last generations of stars to form, 
whereas the stellar metallicity represents an average over 
the entire star formation history of the galaxy. Notably, the 
large scatter in stellar metallicity at fixed gas-phase oxygen 
abundance (solid lines) cannot be solely accounted for by 
the errors in stellar metallicity estimates (indicated by the 
dashed lines). 

To investigate the origin of this scatter, we compute 
the residuals in \og(Z / Zq) with respect to a simple lin- 
ear fit of the median relation (large dots) in Fig. [§] The 
residuals are plotted in Fig. 1101 as a function of various 
galaxy properties, inclu ding the star formation rate (from 
iBrinchmann et al"1l2004l) and an indirect estimate of the gas 
mass fraction inferred from the star formation rate (see eq. 5 



14 A. Gallazzi et al. 



gas ejection and/or accretion may be important factors in 
galaxy chemical evolution. 
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Figure 9. Stellar metallicity estimates (present work) 
plotted against gas- phase oxygen abundance 12+log(0/H) 
jTremonti e t al. 2004) for 7462 high-S/N SDSS-DR2 galaxies 
for which both measures are available. The small dots show 
the median-likelihood estimates of the two parameters for each 
galaxy. The larger points (joined by a solid line) show the median 
stellar metallicity in bins of 12+log(0/H) (each bin containing 
~ 300 galaxies), while the outer solid lines show the corre- 
sponding 16 th and 84 th percentiles. The dashed lines indicate 
the mean 68 percent confidence range in the stellar metallicity 
estimates for the galaxies in e ach bin of 12+log(Q/H). Solar 
metallicity is 8.69 in these units lAllende Prieto et al 120011) . The 
dot-dashed line represents the one-to-one relation, showing that 
stellar metallicity is always lower than gas-phase metallicity. 

of lTremonti et alJl2004l) . In each panel, we indicate the me- 
dian trend of the residuals as a function of the property un- 
der consideration for galaxies in five bins of gas-phase oxygen 
abundance, from 8.70±0.075 (blue) to 9.3±0.075 (red). The 
\og{Z / Zq) residuals show positive correlations with stellar 
mass, surface mass density, concentration parameter and 
D4000 and negative correlations with specific star forma- 
tion rate, gas mass fraction and H0A+H7A. These various 
trends are almost independent of gas-phase oxygen abun- 
dance. 11 They imply that, at fixed 12+log(0/H), the stars 
are most metal-rich in massive, high-concentration galaxies, 
which have exhausted most of their gas and form stars at 
a lower rate than in the past. The existence of a substan- 
tial range of stellar metallicities at fixed gas-phase oxygen 
abundance in Fig. [§] further suggests that the galaxies in 
our sample are not well approximated by 'closed-box' sys- 
tems, for which stellar and interstella r metallicities should 
be tightly related (e.g.. lTinslevlll98d) . This indicates that 



11 We checked that the trends of \og(Z / Zq) residuals with the 
strengths of D4000 and +H7.A are not caused by any bias in 
our estimates of stellar metallicity (see Section l2.4.2l . 



3.3 Age versus metallicity 

Most previous population synthesis studies of the age- 
metallicity relation fo r nearby galaxies have focused on 
early-type galaxies (e . g. |Gonzalez et alll99"3lWorthey|l997l; 
[Bernardi et alj |l99&t [ Terlevich et alJ Il999l : iFerreras et alJ 
Il999t iTraeer et alJ 2000). A few studies have also in- 
clude d the bulges of spira l galaxies (lJabJonk aet_ap 
I l996t iGoudfrooii et all Il999t iTerlevich fe Forbes! 120021 
iProctor fe Sansoml 12002 j. Our determination, with a new 
method, of the ages and metallicities of a large sample of 
nearby galaxies spanning a wide range of star formation ac- 
tivities allows us to re-assess the relation between these two 
physical parameters. 

We separate here high-concentration, bulge-dominated 
early-type galaxies and low-concentration, disk-dominated 
late-type galaxies as described above to understand whether 
different conclusions can be drawn for different types of 
galaxies. We have shown in Fig. |5] that both metallicity 
and age increase with stellar mass. We may therefore ex- 
pect that the distribution of galaxies over metallicity and 
age will change when different masses are considered. This 
mass dependence may also be different for different types of 
galaxies. For this reason, we divide each morphological sub- 
sample into six bins of stellar mass. Figs lf ll and !12l show the 
age-metallicity relation in these different stellar-mass bins 
for late-type and early-type galaxies, respectively. The stel- 
lar mass increases from log(M/M0) ^ 10 in the top-left 
panel to log(A4/MQ) > 11 in the bottom-right panel. These 
relations were obtained by simply coadding the normalized 
joint likelihood distributions of age and metallicity for in- 
dividual galaxies, without further normalization. The grey 
levels indicate the probability associated with each value of 
metallicity and age in a square-root scale, while the contours 
enclose 26, 68, 95 and 99 percent of the total probability. 

Fig. 1111 shows that the ages and metallicities of low- 
concentration, late-type galaxies change significantly as a 
function of stellar mass. Galaxies less massive than 10 10 M© 
have typically low metallicities [Z « 10~ 6 Zq) and young 
ages (0.8 < i, /Gyr < 2.5). At intermediate stellar masses, 
i.e. for 10 10 < M*/M Q < 10 10 5 , the distribution in metal- 
licity becomes broader and shifts gradually to higher metal- 
licities. The overall spread in age is roughly 0.8 dex, and the 
distribution shifts slightly to older ages. Late-type galax- 
ies more massive than 1O 1O ' 5 M0 have predominantly high 
metallicities (Z > Zq) and old ages (3.2 < i r /Gyr < 10). 
The fraction of young metal-poor galaxies nearly vanishes 
at the highest masses. 

Fig. 1121 shows that the ages and metallicities of 
high-concentration, early-type galaxies depend much more 
weakly on stellar mass than those of low-concentration, late- 
type galaxies. The gradual increase of metallicity with stel- 
lar mass is still appreciable for early-type galaxies. Galax- 
ies less massive than 10 10 Mq have metallicities typically in 
the range 1O"°' 4 -1O O1 Z , while galaxies more massive than 
10 11 M© have metallicities typically above solar. The age 
range does not exhibit any significant variation with stellar 
mass, except for a population of young metal-poor galaxies 
in the lowest-mass bin, which does not show up at higher 
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Figure 10. Residuals in log(Z/ZQ) with respect to a simple linear fit of the median relation between stellar metallicity and gas-phase 
oxygen abundance in Fig. plotted against various galaxy properties. From left to right, top raw: star formation rate (in units of 
MQyr -1 ), specific star formation rate (in units of yr — 1 ) and gas mass fraction. Middle raw: r-band light-weighted age (in units of yr), 
D4000 and Hi^+IPm index strengths. Bottom raw: total stellar mass (in units of Mq), surface mass density within the Petrosian r-band 
-R50 radius (in units of Af0kpc — 2 ) and r-band concentration parameter. In each panel, curves of different colors indicate the median 
relations in different bins of gas-phase oxygen abundance, increasing from 8.7 ± 0.075 (blue) to 9.3 ± 0.075 (red). 



masses. At low stellar masses, we also notice a small fraction 
(~ 1 percent) of galaxies clumping at much higher metallici- 
ties and younger ages than the bulk of the early-type galaxy 
sample. Detailed inspection reveals that these galaxies are 
photometric outliers, for which we are uncertain of the ac- 
curacy of our fits. Follow-up observations would be required 
to understand whether the unusual metallicities and ages 
of these galaxies are consistent with a substantial recent 
episode of metal-enriched star formation. 

We note that the typical light- weighted age for the mas- 
sive early-type galaxies in our sample is ~8 Gyr, which is 
somewhat younger than the ages usually quoted for (clus- 
ter) ellipticals. Differences can arise from differences in the 
adopted methods. In particular, the ages usually quoted in 
the literature are derived by comparing the observed dis- 
tributions in spectral indices with simple stellar population 
predictions. What we quote here is a proper light-weighted 
age, calculated by weighting the age of each generation of 
stars by their luminosity along the entire SFH of each model 
galaxy. Nevertheless, the relatively young ages we find could 
reflect an environmental dependence, since the SDSS sam- 
ple is composed predominantly of galaxies in lower density 
regions than those used in most studies of early-type galax- 



ies. Indeed, several studies have reported an age difference 
between early-type ga laxies in clusters and in low-density re- 
gions from ~1.2 Gvr |Be rnardi et al]ll9 98l) u p to 2 — 3 Gyr 
jKuntschner et alj r2002 ; Tho mas et al]H)05). We will dis- 
cuss the properties of the early-type galaxies in our sample 
in more detail in the second paper of this series. 

A comparison of Figs II II and 1121 suggests that the ages 
and metallicities of late- and early-type galaxies differ pri- 
marily at the low-mass end. Low-mass late-type galaxies 
tend to be younger and more metal-poor than their early- 
type counterparts (by about 0.6 dex in both age and metal- 
licity). In contrast, massive late-type galaxies tend to have 
old ages and high metallicities similar to those of their early- 
type counterparts, even though there is a tail of young 
(tr < 10 9 ' 5 yr), metal-poor [Z « 10~°' 5 Z Q ) late-type galax- 
ies with masses greater than lO 11 A/0 in our sample. The 
high metallicities of most massive late-type galaxies might 
reflect the influence of their metal-rich bulges. This may be 
enhanced by the fact that the SDSS spectra sample only the 
inner regions of the galaxies (see Section f3.4l below) . 

Figs II II and 1121 also indicate that the relation between 
stellar metallicity and stellar mass discussed in Section f3. 21 
holds for both late-type and early-type galaxies. In con- 
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trast, the relation between age and stellar mass is much 
more pronounced for late-type galaxies, which cover a sig- 
nificant range in age (~ 9 Gyr). Early- type galaxies do 
not show such a clear trend, at least above 1O 1O ' 5 M0. 
This is in agreement with previous findings of no signifi- 
cant relatio n between the ages of early-type galaxies and 
luminosity (iKun tschner fc Daviesl If 9981 : iKuntschnerl If 99&I 
iTerlevich fc Forbesll2002jr At fainter luminosities, however, 
there are indications of a spread to younger ages (see also 
IWorthev fc Ottavianll99l : lcaldwell fc Rosell99sh . We also 
find that the fraction of young early-type galaxies increases 
at lower stellar masses. This is true in particular below 
lO^ M^, reflecting a possib le contamination by SO galaxies 
(e.g. IPoggianti et al]l200ll) . Finally, we note that for both 
early-type galaxies and high-mass late-type galaxies, the 
probability contours of Figs If f I and If 21 hint at an anticorre- 
lation between age and stellar metallicity, at given mass. 
Other workers have pointed out this an t icorrelation (e.g . 
Worthev et aH If 9951: IColless et al l If 9991: IJorgensenl If 999 



Table 3. Average slope of the age-metallicity degeneracy in dif- 
ferent bins of stellar mass, for the low-concentration and high- 
concentration galaxies of Figs llll and ll2l 



Rakos et alJl200lt IPoggianti et al.ll200ll) . but it is difficult 



to assess its significance here because of the correlated er- 
rors on age and metallicity (see below). 

The dependence of the properties of SDSS galaxies on 



stellar mass has 


seen addressed in previous studies (e.g. 


Shen et alJ 120031 


Brinchmann et al.l 120041: iTremonti et al. 


2004 


Kannaopan 


12004 ISalim et alJ 120051: Ijimenez et al 


2005 


. In particular, 


Kauffmann et al. 


(2003b) point out 



that, at masses around ~ 1O 1O ' 5 M0, galaxies separate into 
two distinct classes of spectral and structural properties. 
IBaldrv et alJ |2004) have recently shown that this transi- 
tion around (2 - 3) x 1O 1O M can be seen directly from the 
colour-magnitude relation. Moreover, this mass range seems 
also to correspond to a shift in gas richness llKannappanl 
l2004f) . It is interesting to note that, around the same value 
of stellar mass, there appears to be also a gradual transition 
in the ages and metallicities of the galaxies. This transition, 
which is particularly evident for late-type galaxies, causes 
the dispersion in age and metallicity to be highest for galax- 
ies of intermediate mass. 

We have not mentioned so far the potential influence 
of the 'age-metallicity degeneracy' on the results of Figs llll 
and If 21 fsee Section 0. This degeneracy, which arises from 
the competing effects of age and metallicity on absorption- 
line strengths, can cause the confusion bet ween old, meta l- 
poor galaxies and young, metal-rich ones IWorthev 1994). 
To address this issue, we first evaluate the slope of the age- 
metallicity degeneracy from the two-dimensional likelihood 
distribution of each individual galaxy in our sample. This 
is achieved by computing the matrix of second moments of 
the likelihood distribution (inside the 68 percent probability 
contour) with respect to the median log Z and log t r values. 
The direction of the eigenvector corresponding to the largest 
eigenvalue of this matrix defines the slope A log Zj A log t r 
of the age-metallicity degeneracy. In each panel of Figs If II 
and !12l the straight line indicates the mean A log Zj A log t r 
for the galaxies falling into the corresponding bin of stellar 
mass. The slopes are listed in Table OP They tend to be 



Stellar mass (Mq) 



A log Z/A log t r 
C>2.i 



log Ah ^ 10 


-1.78 


-1.29 


10 < log Ah 10.3 


-1.41 


-0.85 


10.3 < logM, ^ 10.5 


-1.28 


-0.84 


10.5 < logM, ^ 10.8 


-1.15 


-0.80 


10.8 < log Ah «5 11 


-1.05 


-0.78 


logM, > 11 


-0.95 


-0.75 



steeper for low-concentration galaxies (—1.78 to —0.95) than 
for high-concentration ones (—1.29 to —0.75). 

The orientation of the age-metallicity degeneracy in 
Figs llll and ll2l indicates that it cannot be responsible for the 
trends in age and metallicity as a function of stellar mass 
identified above. However, it can account for the shape of 
the joint age-metallicity probability contours for early-type 
galaxies and the most massive late-type galaxies at high 
metallicities (Z > Zq). 



3.4 Aperture bias 

The SDSS galaxy spectra on which our analysis is based 
are taken using 3"-diameter fibres, which collect only the 
light coming from limited inner regions of the galaxies. We 
can estimate the fraction of the total galaxy light enter- 
ing a fibre from the ratio between the fibre flux and the 
total Petrosian flux. Fig. If 31 shows the distribution of this 
ratio for the 44,254 high-S/N galaxies in our sample. The 
median value is about 30 percent, and the maximum value 
about 60 percent. The dashed and dotted lines show the 
distributions for low-concentration, late-type galaxies and 
high-concentration, early-type galaxies, respectively. As ex- 
pected, the fibre tends to collect a smaller fraction of the 
total light for low-concentration galaxies than for high- 
concentration galaxies. 

Radial gradients in the properties of stellar pop- 
ulations are known to exist for both bulge-dominated 
jHenrv fc Worthevlll999l: ISaglia et alJ 120001: iMehlert et all 
| 2003ft and disk-do minated jBell fc de Jond l200d: 
iMacArthur et al]|2004l) galaxies. Thus, aperture effects are 
a significant concern for our age and metallicity estimates. 
Fortunately, the distribution of the fraction of total galaxy 
light that enters an SDSS fibre does not vary appreciably 
with redshift (except at very low redshift), because larger 
galaxies are also brighter and hence selected out to larger 
distances. 

One way to test for the effects of aperture bias on our 
results is to compare the ages and metallicities of galaxies 
of similar luminosity located at different redshifts. Because 



These values are consistent with the results of lWorthevHl99d) . 
which imply that the change in metallicity needed to compensate 
a change in age at fixed index strength in the spectrum of an old 



stellar population should correspond to values of A log Z/A log t r 
ranging from —1.7 for age-sensitive indices to —0.2 for metal- 
sensitive indices. 
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Figure 11. Two-dimensional distribution of stellar mctallicity and age in 6 bins of stellar mass (as indicated) for the subsample of 5616 
late-type (C 2.4) SDSS-DR2 galaxies with median S/N per pixel greater than 20. Each distribution was obtained by coadding the 
normalized joint likelihood distributions of age and metallicity for individual galaxies. The grey levels indicate the probability associated 
with each value of metallicity and age in a square-root scale, while the contours enclose 26, 68, 95 and 99 percent of the total probability. 
The straight line (which intersects the peak of the probability distribution) indicates the mean slope of the age-metallicity degeneracy 
for the galaxies falling into each stellar-mass bin. 
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Figure 12. Same as Fig. 1111 but for the subsample of 26,003 early-type (C J; 2.8) SDSS-DR2 galaxies with median S/N per pixel greater 
than 20. 
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Figure 13. Distribution of the fraction of the total r-band flux 
collected by the SDSS fibre, computed as the ratio between the 
fibre flux and the Petrosian flux, for 44,254 high-S/N galaxies in 
the SDSS DR2. The vertical bar indicates the median fraction for 
this sample. The dashed and dotted lines show the distributions 
of low-concentration, late-type galaxies and high-concentration, 
early-type galaxies, respectively. 



of the presence of both a faint and a bright magnitude cut- 
off in the SDSS survey, this comparison is possible only for 
galaxies in a limited absolute- luminosity range, if we want to 
span a reasonably large redshift range. We consider galax- 
ies with r-band absolute magnitudes in the narrow range 
—21 ^ r —22, which are observed at redshifts between 
0.02 and 0.12. At fixed absolute magnitude within this range, 
both low- and high-concentration galaxies show a decrease 
(by ~ 0.15 dex) in typical metallicity from redshift 0.02 to 
0.12. High-concentration galaxies do not exhibit any signif- 
icant trend in age, but the typical age of low-concentration 
galaxies decreases by about 0.2 dex over the same redshift 
range. 

Another way to examine the potential effects of aper- 
ture bias on our results is to compare the age and metallicity 
distributions of galaxies with similar properties as a func- 
tion of z/z m „, where z max is the minimum of the redshift 
at which the galaxy reaches the survey faint-magnitude limit 
and that at which its median S/N per pixel drops below 20 
(we assume Poissonian noise, i.e. that the median S/N per 
pixel scales as the square root of the apparent fibre luminos- 
ity). With this approach, we can avoid restricting our test 
to galaxies in a narrow magnitude range. Also, since we are 
interested in possible effects of aperture bias on trends with 
galaxy mass, we divide our sample into the same 6 bins of 
stellar mass as in Figs II II and 1121 above. 

Fig. IMI shows the dependence on z/z max of the median 
(thick line) and the 16* and 84 th percentiles (thin lines) of 
the distributions in stellar metallicity (left-hand plot) and 
age (right-hand plot) of the high-S/N galaxies in our sam- 
ple. The solid lines correspond to the subsample of 26,003 
early-type galaxies, while the dot-dashed lines correspond 



to the subsample of 5616 late-type galaxies. For early-type 
galaxies of different masses, the largest change in median 
metallicity from one edge of the survey to the other is about 
0.2 dex. This is comparable to the error associated with 
the metallicity estimates of most galaxies and smaller than 
the overall metallicity change with galaxy mass (see Fig. [SJ . 
This variation is also consistent with the metallicity gradient 
of —0.2 dex per decade in radius r eported by several stud- 
ies of bulge-dom inated galaxies fsee lHenrv fc Worthevll999l : 
IWu et al.l20"04Tl . We note that the largest variation arises for 
galaxies with intermediate stellar masses, which often have 
both a prominent bulge and a disk. For disk-dominated, late- 
type galaxies, the trends in median metallicity with z/z max 
are weaker than for early-type galaxies, although the distri- 
butions are broader. 

The median age of early-type galaxies does not show 
any significant trend with z/z max in any mass range in 
Fig. 1141 This is in agreement with ea rlier findings of neg- 
ligible age grad ients in such galaxies (IMehlert et al.ll2003l : 
IWu et al-lEool . Late-type galaxies of low and intermediate 
mass show a similar behaviour. Only the most massive late- 
type galaxies show a significant decrease (by about 0.2 dex) 
in median age as a function of z/z max , as expected if age 
gradients in spiral gala xies are stronger for larger galaxies 
jMac Arthur et aljEooi'l . 

The occurrence of negative radial age and metallicity 
gradients in different types of galaxies and the fact that the 
SDSS spectra sample no more than 50-60 percent of the to- 
tal light of a galaxy are likely to cause overestimates of the 
stellar metallicities and ages of some galaxies in our sample. 
However, the weak trends in log Z and log t r as a function of 
z I Zmax in Fig. 1141 suggest that aperture biases are unlikely 
to have a significant effect on the trends we have found pre- 
viously for metallicity and age as a function of stellar mass. 



4 SUMMARY AND CONCLUSIONS 

We have used a new approach to derive estimates of light- 
weighted stellar metallicity, age and stellar mass from the 
optical spectra of a sample of ~ 200,000 nearby galax- 
ies (0.005 < z < 0.22) drawn from the SDSS DR2. Our 
method relies on the comparison of the galaxy spectra to 
a large library of model spectra at medium-high resolution, 
based on Monte Carlo star formation histories spanning the 
full physically plausible range. Extending earlier work by 
iKauffmann et alJ J2003al) . we have adopted a Bayesian ap- 
proach to derive the a posteriori likelihood distribution of 
each physical parameter by computing the goodness of fit of 
the observed spectrum for all the models in the library. In 
practice, we compute only the strengths of a set of carefully 
selected spectral absorption features. 

Our analysis shows that stellar absorption features with 
different sensitivities to age and metallicity must be fitted 
simultaneously in galaxy spectra to obtain good constraints 
on both age and metallicity. We focus on several newly cal- 
ibrated Lick indices and the 4000 A break (see below). The 
resolution of our models is higher than that used in most 
previous studies and is well matched to the resolution of the 
SDSS spectra. This allows us to measure the indices in the 
same way in model and galaxy spectra, avoiding any loss of 
information from the observed spectra. We emphasize that, 
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Figure 14. Dependence on z/2 max of the median (thick line) and the 16 th and 84 th percentiles (thin lines) of the distributions in stellar 
mctallicity (left-hand plot) and age (right-hand plot) of SDSS-DR2 galaxies with median S/N per pixel greater than 20. The quantity 
z mH is the minimum of the redshift at which the galaxy reaches the survey magnitude limit and that at which its median S/N per 
pixel drops below 20. The solid lines correspond to the subsample of 26,003 early-type (C J; 2.8) galaxies, while the dot-dashed lines 
correspond to the subsample of 5616 late-type (C $C 2.4) galaxies. 



while most previous studies were restricted to the analy- 
sis of old stellar populations, our sample includes galaxies 
in the full range of star formation activities. For actively 
star-forming galaxies, the higher resolution of the models is 
crucial to separate strong emission lines from the underlying 
stellar absorption. 

Using this approach, we are able to constrain stellar 
metallicity and light-weighted age to within ±0.15 dex for 
the majority of the galaxies in the sample. Our stellar-mass 
estimates, determined within less than ±0.1 dex for the 
whole sample, are in good agreement with those derived by 
Kau ffmann et all (l2003al) . While the uncertainties on the 
stellar mass are almost independent of galaxy type, we find 
that, as expected, the smallest errors on age and especially 
metallicity are obtained for galaxies with the strongest ab- 
sorption lines. 

The estimates we derive of galaxy parameters are af- 
fected by limitations of both model and observed spectra. 
For example, the models do not include the effects of vari- 
ations in heavy-element abundance ratios, while systematic 
deviations from the solar abu ndance ratios are kn own to 
arise in external galaxies (e.g. IWorthev et alJll993) . When 
estimating metallicity, we therefore consider only those spec- 
tral features that have been shown by previo us studies to 
depend negligibly on q/F e abundance ratio iTrager et alJ 
l2000l: iThomas et alll2003bl iTantalo fc Chiosill2004l) . The in- 
clusion of higher-order Balmer lines, which may depend on 
eleme nt abundance ratio at high metallicities jThomas et alJ 
120041) . could lead to an overestimate of the metallicity and 
an underestimate of the age by ~0.05 dex for galaxies with 
non-solar abundance ratios in our sample. Another system- 
atic uncertainty in our metallicity and age estimates comes 
from the discrepancy between integrated colours and spec- 
tral index strengths for supra-solar metallicity models (the 



model spectral indices correspond to lower metallicities than 
the colours). Therefore, we expect metallicities inferred us- 
ing these models to be biased toward high values at high 
metallicities. This problem mainly affects massive, early- 
type galaxies with strong absorption lines. 

We have shown that, in addition to these model un- 
certainties, the signal-to-noise ratio in a galaxy spectrum 
strongly affects our estimates of age and especially metallic- 
ity. We choose to rely primarily on results for galaxies with 
median S/N per pixel greater than 20. This cut reduces the 
original sample by about 75 percent and biases our analysis 
toward high surface brightness, high concentration galaxies. 
Our main results appear unaffected by this cut in S/N; in- 
cluding lower-S/N galaxies mainly increases the scatter in 
the relations we find, because of the larger associated uncer- 
tainties on metallicity. 

Another observational limitation is the small aperture 
sampled by the SDSS spectra. The fraction of light collected 
by an SDSS fibre depends on the apparent size of the target 
galaxy, and this fraction is less than 30 percent for most 
galaxies in our sample. The stellar mass we derive from 
the spectrum has been scaled to the total luminosity of the 
galaxy, and hence, it represents the total stellar mass of that 
galaxy. However, the age and the metallicity derived from 
the fibre spectrum cannot be easily extrapolated to total 
values without information about radial gradients. We have 
shown that these two parameters vary weakly as a func- 
tion of the fraction of light entering the fibre. This reassures 
us that aperture bias does not have a major effect on the 
relations between metallicity, age and stellar mass that we 
derive. 

We have explored how galaxies are distributed in 
metallicity, age and stellar mass in the diagram defined 
by 4000-A break and Balmer-line strength. As shown by 
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Kauffmann ct al ] J2003al) the principal sequence in this 
plane reflects an increase in present to past-averaged 
star formation rate w ith decreasing D4000 (see also 
iBrinchmann et al"ll2004f) . Galaxies with particularly strong 
Balmer absorption at given D4000 have undergone recent 
starbursts. By including metallicity dependence in the mod- 
elling, we have shown that: 

• In this diagram, the average stellar mass increases along 
the sequence populated by galaxies in our sample, from less 
than 1O 1O M0 to more than lO n Af0 from end to end of the 
sequence. The transition at intermediate D4000 is smooth 
between the two regimes. These resu lts are in agreement 
with those of iKauffmann et"al 

• Both stellar metallicity and light-weighted age also in- 
crease with D4000. Galaxies with particularly strong Balmer 
absorption at intermediate 4000-A breaks (D4000 ~ 1.5) 
are on average younger and more metal-rich than the bulk 
of the sequence. This supports the idea that these galaxies 
may have recently undergone a burst of metal-enriched star 
formation. 

In addition, we have constructed the conditional dis- 
tributions of metallicity and age estimates as a function of 
stellar mass. This analysis shows that: 

• Both stellar metallicity and light-weighted age increase 
with stellar mass, the increase being rapid at intermediate 
masses. At masses above ~ 3 x 1O 1O M0, a gradual flattening 
occurs in both relations. This st ellar mass corresponds to 
the transition mass identified bv IKauffmann et alJ (|2003b) 
in plots of D4000, surface mass density and concentration 
against stellar mass. 

• Despite the above clear relations, metallicity and age 
are not uniquely determined by mass; there is an intrinsic 
scatter in both relations, which is largest at stellar masses 
around 10 10 Mq. This scatter persists even when considering 
low- and high-concentration galaxies separately. 

• The relation we find between stel lar metallicity and stel- 
lar mass is similar to that found bv lTremonti et al] J2004T) 
between gas-phase oxygen abundance and stellar mass. We 
confirm that higher stellar metallicities are indeed associated 
with higher gas-phase metallicities. However, the relation 
between stellar and gas-phase metallicities has a substantial 
scatter, which can only partly be attributed to the uncer- 
tainties in stellar-metallicity estimates. This suggests that a 
simple 'closed-box' scenario is not sufficient to explain the 
relation, and that a variety of gas accretion/ejection histo- 
ries may be required. 

We have also explored relationships between metallic- 
ity, age and stellar mass for separate subsamples of high- 
concentration, early-type and low-concentration, late-type 
galaxies. The distribution in age-metallicity space in bins of 
stellar mass shows that (Figs fTTl and 1121 : 

• Relations between stellar metallicity, age and stellar 
mass hold for both early- and late-type galaxies. However, 
the age range of the sample is less mass-dependent for early- 
type than for late-type galaxies. 

• At masses below 1O 1O M0, late- type galaxies are younger 
and more metal-poor than early-type galaxies. At masses 
above 1O 11 M0 the metallicities and ages of late-type galaxies 
are similar to those of early-type galaxies. 



The above results indicate that young, metal-poor stel- 
lar populations are found predominantly in low-mass galax- 
ies. This is consistent with a 'downsizing' scenario, in which 
the mass and luminosity of the galaxies undergoing active 
star formation become progressiv ely lower as the Universe 
becomes older llCowie et al]ll996h . Taken at face value, this 
would imply that the low metallicities of low-mass galax- 
ies just reflect the fact that these galaxies have had less 
time to produce metals. This naiv e conc lusion is not consis- 
tent with that of iTremonti et all <l2004f) . who showed that 
the effective yields (from indirect gas-mass fraction measure- 
ments) of low-mass, star-forming galaxies in the SDSS DR2 
are too low for these galax ies to be interp r eted as young 
closed-box systems. This led ITremonti et all <l2004f) to favor 
galatic winds as the most likely origin of the relation between 
gas-phase metallicity and stellar mass. This interpretation is 
supported here by the fact that galaxies with no Ha emission 
in our sample (which have presumably completed their star 
formation) also show a trend of decreasing stellar metallic- 
ity with decreasing stellar mass for M, < 2 x 10 10 Af Q (the 
metallicity of closed-box galaxies should tend to the yield 
near gas exhaustion, independently of mass). The large scat- 
ter present in all relations between stellar metallicity, age 
and stellar mass (even when separated into late-type and 
early-type galaxies) further indicates that stellar mass is not 
the unique parameter determining the star formation history 
and hence the physical parameters of present-day galaxies. 
Additional variations in star formation and enrichment his- 
tory (e.g. driven by gas infall or outflow) are required to 
explain these parameters. 

An advantage of the large samples available from the 
SDSS is that they allow us to study not only the mean re- 
lations between physical parameters of galaxies, but also 
the dependence of these parameters on various galaxy prop- 
erties. Our results are particularly useful for studying the 
physical origin of well-known observed relations, such as for 
example the colour-magnitude relation and the relation be- 
tween Mg2-index strength and velocity dispersion for early- 
type galaxies. Our analysis also enables us to estimate the 
total metal content of the local Universe and the distribution 
of these metals over different galaxy types. We will address 
these questions in forthcoming papers. 

A complete census of the physical parameters of galax- 
ies out to z ~ 0.1 is essential for constraining models of the 
star formation and chemical enrichment histories of galaxies. 
The application of our method to large samples of galaxies 
at higher redshifts will enable us to study how the distri- 
bution of the physical parameters of galaxies evolves with 
lookback time. O ngoing deep redshi ft surveys, such as VVDS 
<Le Fevre et al.l2004l) and DEEP2 JPavis et all2003h in th e 
optical and GALEX in the ultraviolet jMartin et afll2004h . 
are already assembling large samples of high-redshift galax- 
ies. Quantitative comparison with the SDSS samples should 
substantially deepen our understanding of galaxy evolution. 
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